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ABSTRACT 

Measurements of small-scale turbulent fluctuations in the solar wind find a non-zero right-handed 
magnetic helicity. This has been interpreted as evidence for ion cyclotron damping. However, theoret- 
ical and empirical evidence suggests that the majority of the energy in solar wind turbulence resides in 
low frequency anisotropic kinetic Alfven wave fluctuations that are not subject to ion cyclotron damp- 
ing. We demonstrate that a dissipation range comprised of kinetic Alfven waves also produces a net 
right-handed fluctuating magnetic helicity signature consistent with observations. Thus, the observed 
magnetic helicity signature does not necessarily imply that ion cyclotron damping is energetically 
important in the solar wind. 
Subject headings: turbulence — solar wind 



1. INTRODUCTION 

The identification of the physical mechanisms respon- 
sible for the dissipation of turbulence in the solar wind, 
and for the resulting heating of the solar wind plasma, re- 
mains an important and unsolved problem of hcliosphcric 
physics. An important clue to this problem is the ob- 
served non-zero fluctuating magnetic helicity signature 
at scales corresponding to the dissipation range of solar 
w ind turbulence. 

iMatthaeus et alJ (|1982f ) first proposed the "fluctuat- 
ing" magnetic helicity as a diagnostic of solar wind tur- 
bulence, defining the "reduced fluctuating" magnetic he- 
licity spectrum derivable from observational data (see 
$3] below) . A subsequent study, corresponding to scales 
within the inertial range, found values that fluctuated 
randomly in sign, and suggested an interpretation that 
"a substantial degree of helicity or circular polarization 
exists throughout the wavenumber spectrum, but the 
sense of polarization or handed ness alternates randomly" 
(jMatthaeus fc Gold stein 1982ft . Based on a study of the 
fluctuating magnetic helicity of solu tions to the linear 
Vlasov-Maxwell dispersion relation, iGarvl (|1986h sug- 
gested instead that, at inertial range scales, all eigen- 
modes have a very small intrinsic normalized fluctuating 
magnetic helicity, eliminating the need to invoke an en- 
semble of waves with both left- and right-handed helicity 
to explain the observations. 

Subsequent higher time resolution measurements, cor- 
responding to scales in the dissipation range, exhibited a 
non-zero net reduced fluctuating magnetic helicity signa- 
ture, with the sign apparent ly correlated with the direc- 
tion of the magnetic sector ((Goldstein et al.|[l994f ). As- 
suming dominantly anti-sunward propagating waves, the 
study concluded that these fluctuations had right-handed 
helicity. The proposed interpretation was that left-hand 
polarized Alfven/ion cyclotron waves were preferentially 
damped by cyclotron resonance with the ions, leaving un- 
damped right-hand polarized fast/ whistler waves as the 



dominant wave mode in the dissipation range, producing 
the measured net reduced fluctuating magnetic helicity. 
We refer to this as the cyclotron damping interpretation. 

A subsequent analysis of more solar wind inter- 
vals confirmed these fin dings for the d i ssipatio n range 
(jLeamon et alJ fl998bf ). iLeamon et al.l (|1998af ) argued 
that a comparison of the normalized cross-helicity in the 
inertial range (as a proxy for the dominant wave prop- 
agation direction in the dissipation range) to the mea- 
sured normalized reduced fluctuating magnetic helicity 
provides evidence for the importance of ion cyclotron 
damping, which would selectively remove the left-hand 
polarized waves from the turbulence; using a simple rate 
balance calculation, they concluded that the ratio of 
damping by cyclotron resonant to non-cyclotron reso- 
nant dissipation mechanisms was of order unity. A re- 
cent study performing the same analysis on a much larger 
data set concurred with this conclusion ([Hamilton et al.l 
l200l . 

In this Letter, we demonstrate that a dissipation range 
comprised of kinetic Alfven waves produces a reduced 
fluctuating magnetic helicity signature consistent with 
observations. A dissipation range of this nature re- 
sults from an anisotropic cascade to high perpendicu- 
lar wavenumber with k± 3> kn ; such a cascade is con- 
sistent with existing theories for low-f r equency plasma 
turbulence (iGoldreich & Sridharl Il995|: iBoldvrevi 120061 : 
iHowes et al.ll2008bt iSchekochihin et al.ll2009fl, numerical 
simulations ( Cho &; Vishniad 120001 : Howes et al.l 12 008a'). 



and observati ons in the solar wind ( Horburv et alJ 20081 : 
Podesta 2009) . Our results imply that no conclusions can 
be drawn about the importance of ion cyclotron damping 
in the solar wind based on the observed magnetic helicity 
signature alone. 

2. FLUCTUATING MAGNETIC HELICITY 

The magnetic helicity is defined as the integral over the 
plasma volume H m = J d 3 rA ■ B, where A is the vector 
potential which defines the magnetic field via B = V x A. 
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This integral is an invariant of ideal Magnetohydrody- 
namics (MHD) in the absence of a mean magnet ic field 
(|Woltier)ll95"8aH bT). iMatthaeus fc Goldstein! (|1982f ) chose 
to set aside the complications associated with the pres- 
ence of a mean magnetic field, defining the fluctuating 
magnetic helicity by H' m — J drSA ■ (5B, where the fluc- 
tuating quantities denoted by S do not include contribu- 
tions from the mean field. 
Modeling the turbulent magnetic field 1 by 

B(r,i) =B z + ^B(k)e J ( k - r -^ (1) 

k 

in a periodic cube of plasma with volume L 3 , we ob- 
tain H' m = L 3 H' m ()s), where the fluctuating mag- 
netic helicity density for each wave vector k is defined 
by i4(k) = A(k) • B*(k). Here B(-k) = B*(k) and 
u(— k) = —uj*(k) are reality conditions and B*(k) is the 
complex conjugate of the Fourier coefficient. Specifying 
the Coulomb gauge V • A = 0, we obtain 

H^)=i BvB \ B;B '=i B ' B \ B:B ' 

= i B x B y -B x By ^ (2) 

k z 

where the components Bj(\i) arise from the eigenfunc- 
tions of the linear wave mode. It is easily shown that 
this result is invariant to rotation of the wave vector k, 
along with its corresponding linear eigenfunction, about 
the direction of the mean magnetic field. The normalized 
fluctuating magnetic helicity density is defined by 

a m (k) = fcJ4(k)/|B(k)| 2 , (3) 

where k = |k|. This normalized measure has values 
within the range [— 1,+1], where negative values de- 
note left-handed helicity and positive values denote right- 
handed helicity. 

We numerically calculate cr m (k) over the k±— ku plane 
for the eigenmo des of the linear Vlasov-Maxwell dis- 
persion relation (|Stixl I1992T ) for a proton and electron 
plasma with an isotropic Maxwellian equilibrium dis- 
tribution function for each s pecies and no drift veloc- 
ities (see iHowes et al.l I2006L for a description of the 
code). The dispersion relation depends on five param- 
eters lo = LO VM {k^Pi,k\\p t , Pi,Ti/T e ,v t hJ c), for ion Lar- 
mor radius pi, ion plasma beta ft, ion to electron temper- 
ature ratio Tj/T e , and ion thermal velocity to the speed 
of light v t h z I c 

We specify plasma parameters characteristic of the so- 
lar wind at 1 AU: ft = 1, T,/T e = 1, and v th Jc = 1(T 4 . 
Figure [T] is a contour plot of <r TO (k) obtained by solving 
for the Alfven wave root over the k±—k\\ plane, then us- 
ing the complex eigenfunctions to determine a m (k) . The 
MHD regime corresponds to the lower left corner of the 
plot, knpi <C 1 and k±pi <C 1; here, the Alfven wave with 
fcj. ~ fcii is linearly polarized with o~ m ~ 0. As one moves 
up vertically on the plot to the regime fcii 3> fcj_, the so- 
lution becomes left-handed with values of <r m — > — 1 . In 

1 We assume that turbulent fluctuations are reasonably modeled 
as a collection of linear wave modes. Nonlinear interactions, ne- 
glected here, will serve to replenish energy lost from wave modes, 
so we neglect the linear wave damping and take only the real fre- 
quency. 
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Fig. 1. — Normalized fluctuating magnetic helicity density c m (k) 
[eq. [3] for linear Alfven waves over the kj_-k^ plane. The MHD 
Alfven wave (MHD Alfven) , ion cyclotron wave (IC W) , and kinetic 
Alfven wave (KAW) regimes are identified. Plasma parameters are 
representative of the near-Earth solar wind. 

this regime of nearly parallel wave vectors, the solution 
represents Alfven waves in the limit knpi -C and ion 
cyclotron waves in the limit knpi ^ VWi- The linear wave 
mode becomes strongly damped via the ion cyclotron 
resonance at a value of kupt > \f]Ti (|Garv fc Borovskvl 
2004). This is precisely the behavior supporting the cy- 
clotron damping interpretation of the measured magnetic 
helicity in the solar wind. 

But the Alfven wave solution does not always produce 
left-handed magnetic helicity. If one moves instead from 
the MHD regime horizontally to the right, the solution 
becomes right-handed wit h a m — > +1 a s k±pi —> 1, a be- 
havior previously found bv lGarvl (|1986f ). In this regime of 
nearly perpendicular wave vectors with k± 3> fen , the so- 
lution represents Alfven waves in the limit k±pi <C 1 and 
kinetic Alfven waves in the limit k±pi > 1. Thus, if the 
dissipation range is comprised of kinetic Alfven waves, 
as suggested by theories for critically balanced , low- 
frequency plasma t urbulence ([Schekochi hin et al.l 120091 : 
IHowes et al.ll20 08a). one would expect to observe a pos- 
itive normalized fluctuating magnetic helicity signature 
in that regime. 

3. REDUCED FLUCTUATING MAGNETIC HELICITY 

Unfortunately, due to the limitations of single-point 
satellite measurements, equations ([2]) and ([3]) cannot be 
used directly to calculate the fluctuating magnetic he- 
licity from observations; approximations must be intro- 
duced to define a related measurable quantity. In this 
section, we calculate the reduc ed fluctuating m a gnetic 
helicity density, as defined by IMatthaeus et al.l ()1982D 
and used by subsequent authors, for the magnetic field 
defined by equation ([1]) , but without assuming the Taylor 
hypothesis. 

The two-point, two-time magnetic field correlation 
function is 

(r, t) = {6Bi(x, t)5Bj (x + r, r + t)) , (4) 

where the angle brackets specify an ensemble average, 
defined here by (a(r, t)} = L~ 3 J c? 3 xa(x, r, t). We find 

i? y -(r,t) = ^S*(k) J B J (k) e ^ k — *) (5) 

k 
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where the reality conditions ensure that this quantity is 
real. 

We choose to sample this correlation function at a mov- 
ing probe with position given by r = — vt; this corre- 
sponds to satellite measurements of the solar wind, where 
the probe is stationary and the solar wind is stream- 
ing past the probe at velocity v. Thus, we may de- 
termine the reduced magnetic field correlation function, 
i2y(£) = iZjj (r, t) | r= _ vt , obtaining the form 



RUt) 



E 

k 



S*(k)^(k)e- 



-j(k-v+o>)t 



(6) 



The reduced frequency spectrum, defined by Sy(o/) = 
(l/2vr) / dt'R r l0 (t')e w,t ' , is then given by 

sr.( w ') = ]T B;(k)B 3 -(k) v - (k • v + (?) 



This demonstrates that the frequency to' of the fluctua- 
tions sampled by the moving probe is the Doppler shifted 
frequency to' = k • v + to. N ote that adopting the Tay- 
lor hypothesis (|Tavlorl [l938f ). as often done in studies 
of solar wind turbulence, corresponds to dropping u> in 
equation 0. 

The reduced fluctuating magnetic helicity density is de- 
fined by 



H-(h) = 2Im[S r 23 (h)]/h 



(8) 



where the effective wavenumber is calculated from the 
measured frequency using k\ = u)'/v, assuming the 
Taylor hypothesis is satisfie d (|Matthaeus et alj Il982t 
iMatthaeus fc Goldstein! IT982I ), and we have chosen an 
orthonormal basis with direction 1 along the direction 
of sampling v = v/|v| and directions 2 and 3 in the 
plane perpendicular to v. The normalized reduced fluc- 
tuating magnetic helicity density is given by a^ n (ki) = 
fciiJ^(fci)/|B(fci)| 2 , where |B(fci)| 2 is the trace power. 

The relation between the reduced fluctuating magnetic 
helicity density H^(ki) and the fluctuating magnetic he- 
licity density H' m (k) can be seen by writing the spectrum 
in terms of the Doppler-shifted frequency lo' instead of 

fci, fl£(a/) = 2Im[S , 2 3 (w')]/(^' /»)• Using equation ft) 
and 2Im[a*fo] = i(ab* — a*b), the reduced fluctuating mag- 
netic helicity density can be written as 



k 



i[B 2 (k)SJ(k)-B5(k)B 3 (k)] 



u'/v 



xS[u>' - (k-v + w)] 



(9) 



Equation (J5J), the experimentally accessible quantity, is 
in terms of the magnetic field measurements in a frame 
defined by the solar wind velocity v. To write this in 
terms of the theoretically calculable H' m {k) (eq. [2]), we 
express the magnetic field components B2 and B3 in the 
x, y, z coordinate system. To do so, define the probe ve- 
locity in spherical coordinates about the direction of the 
mean magnetic field: v = v sin 9 cos + v sin 9 sin <py + 
v cos 9z. The orthonormal basis specified with respect to 



v can be written as 

ei = v = sin 9 cos 4>x + sin 9 sin (j>y + cos 9z 

62 = z x v/|z x v| = — sin 0x + cos 4>y 

£3 = ei x £>2 = — cos0cos</>x — cos#sin</>y + sin6*z. 

(10) 

Finally, we exploit the fact that the solutions of the 
Vlasov-Maxwell dispersion relation depend only on the 
perpendicular and parallel components of the wave vec- 
tor k±_ and kii with respect to the mean magnetic field, 
and not on the angle about the field; thus the eigenfunc- 
tion for a wave vector k = fcj_x+fc||Z can be rotated by an 
angle a about the mean magnetic field to yield the solu- 
tion for any wave vector k' = k± cos ax. + k± sin ay + fc||Z. 
Using the above, the reduced fluctuating magnetic helic- 
ity density H^co') in equation © becomes 



k± sin 9 cos a + k\\ cos 9 



k± sin 9 cos a + fcii cos 9 + ui/v 

k 11 ' 

x5[w'- (k'- v + w)], (11) 

where we have specified the azimuthal angle of the probe 
velocity 4> — without loss of generality. It is clear from 
this equation that all possible wave vectors k' that give 
the same Doppler shifted frequency lo' will contribute 
to the sum for the reduced fluctuating magnetic helicity 
density at the frequency lo' . 

4. DISCUSSION 

Predicting the values of HJ,(uj') for solar wind tur- 
bulence based on equation Ijlip requires understanding 
three issues: the scaling of the magnetic fluctuation spec- 
trum with wavenumber, the imbalance of Alfven wave 
energy fluxes in opposite directions along the mean mag- 
netic field, and the variation of the angle 9 between the 
solar wind velocity v and the mean magnetic field. 

The 1-D magnetic energy spectrum in the solar wind 

typically scales as k 1 5 ^ 3 in the inertial ran ge and k^ in 
the dissipation range, where —2 < p < —4 (|Smith et alj 
12006( 1 and the effective wavenumber is ki = uj'/v. It is 
clear from equation (|11| that, when the plasma frame fre- 
quency lo is negligible, the Doppler-shifted observed fre- 
quency always results in an effective wavenumber k\ < k, 
with equality occurring only when the velocity v is 
aligned with the wave vector k. We assume that, for 
homogeneous turbulence at the dissipation range scales, 
turbulent energy at fixed k± and feii is uniformly spread 
over wave vectors with all possible angles a about the 
mean magnetic field. Because the fluctuation amplitude 
deceases for larger effective wavenumbers, the contribu- 
tion to H^iuj 1 ) is maximum at angle a = 0; for angles a 
yielding a Doppler shift to lower effective wavenumbers 
ki < {k\ + fcjj) 1 / 2 , the higher amplitude fluctuations at 
those lower wavenumbers will contribute more strongly 
to H^uj'). An accurate calculation of the magnetic he- 
licity signature based on equation pip must take into 
account the scaling of the magnetic energy spectrum. 

To compare to (fci ) derived from observ ations (for 
example, see Figure 1 of lLeamon et al.l ()1998bl )). we con- 
struct the normalized quantity 



£ k [|B(k)|7fc]Sk-(k'-v + w )] ■ 



(12) 
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Fig. 2. — Normalized reduced fluctuating magnetic hclicity 
^m(^l) vs - effective wavenumber k\ due to a turbulent spectrum 
of kinetic Alfven waves with 9 = 60° . The solid line corresponds to 
the model 1-D energy spectrum while the dashed line corresponds 
to a fc~ 1 spectrum. 

In evaluating equation (|12p , we assume a model 1-D en- 
ergy spectrum 2 that scales as fe -5 / 3 for kpi <C 1 and 
fc -7 / 3 for kpi ^> 1, consist ent with theories for crit- 
ically balanced turbulence dGoldreich fc Sridharl 119951 : 
iHowes et all I2008ah ISchekochihin et alJ I2009T > and solar 
wind observations (Smith ct al. 2006). In Figure [21 we 
plot u1 m {kx) vs. effective wavenumber k\ — ui'/v for a 
turbulent spectrum filling the MHD Alfven and kinetic 
Alfven wave regimes (fej_ > k\\ and kupi < 1) for /3j = 1, 
T,/T e = 1, vthjc = 10~ 4 , 9 = 60°, and v/v A = 10. The 
contributions to ^(fci) for all angles a of each wave 
vector are collected in 120 logarithmically spaced bins in 
Doppler-shifted frequency. The results are rather insen- 
sitive to the scaling of the 1-D magnetic energy spectrum 
over the range from fc -1 to fc~ 4 . The solid line in Fig- 
ure [2] corresponds to the model spectrum assumed above, 
while the dashed line corresponds to a fc -1 energy spec- 
trum. Figure demonstrates that turbulence consisting 
of Alfven and kinetic Alfven waves produces a positive 
(right-handed) magnetic helicity signature in the dissi- 
pation range at k\pi > 1. 

The analysis presented in Figure[2]considers only waves 
with fen > 0, so all of the waves in the summation in 
equation (fTTj) are traveling in the same direction. If 
there were an equal Alfven wave energy flux in the oppo- 
site direction — a case of balanced energy fluxes, or zero 
cross helicity — the net cr£j(fci) would be zero due to the 
odd symmetry of H' m (\t) in fen. It is often observed, at 
scales corresponding to the inertial range, that the energy 
flux in the anti-sunward direction dominates, leading t o 
a large normalized cross helicity (|Leamon etanil998aj) . 
If this imbalance of energy fluxes persists to the smaller 
scales associated with the dissipation range, a non-zero 
value of CT^(fci) is expecte d. However, the ories of un- 
balanced MHD turbulence ( Chandranl 120081 . and refer- 
ences therein) predict that the turbulence is "pinned" to 
equal values of the oppositely directed energy fluxes at 
the dissipation scale. This implies that, at sufficiently 
high wavenumber fci, the value of crLffci) should asymp- 
tote to zero. Thus, (^(fci) in Figure [2] would likely drop 

2 On 150 X 150 logarithmic gridpoints over fcxPi,fc||Pi 6 
[10~ 3 , 10 2 ], the model weights B 2 as a function of fc = (fc 2 + k?,) 1/2 



to zero more rapidly than shown, leaving a smaller pos- 
itive net a L, (fci ) around fci pj ~ 1, consistent with ob- 
servations (Go dsteinet al.l Il994t iLeamon et al.l Il998bt 
Hamilton et al.ll2008l ). Wedefer a detailed calculation of 
the effects of imbalance to a future paper. 

The angle between Bo and v is likely to vary during 
a measurement; this angle does not typically sample its 
full range < 9 < n, but has some distribution about the 
Parker spiral value. Calculations of o-^(fci) over < 9 < 
7r yield results that are qualitatively similar to Figure [2j 
so this averaging will not significantly change our results. 

Taken together, we have demonstrated that a solar 
wind dissipation range comprised of kinetic Alfven waves 
produces a magnetic helicity signature consistent with 
observations, as presented in Figure [2] The under- 
lying assumption of the cyclotron damping interpreta- 
tion of magnetic helicity measurements, an interpretation 
that d ominates the solar win d literature ([G oldstei n et ail 
119941: ILeamon et all U998aHbl : lHamilton et al.l I2008D . "is 
the slab model, k = fcyz and k± = 0, i.e., purely par- 
allel wave vectors. As shown in Figure [IJ only in the 
limit fen ^> k± does the Alfven wave root generate a left- 
handed helicity a m — > — 1 as knpi — > \f$i\ in the same 
limit, the fast/ whistler root generates a right-handed 
helicity a m — > +1 in a quan titatively similar manner 
(see Figure 9 of iGarvl (|1986l )). Strong ion cyclotron 
damping of the Alfven/ion cyclotron waves as knpi — ► 1 
(jGarv fc Borovskvl 120041 ) would leave a remaining spec- 
trum of right-handed fast/ whistler waves, as proposed by 
cyclotron damping interpretation. However, only if the 
majority of the turbulent fluctuations have feii > k± is the 
slab limit applicable, and only if significant energy resides 
in slab-like fluctuations are the conclusions drawn about 
the importance of cyclotron damping valid. There is, on 
the other hand, strong theoretical and empirical support 
for the hypothesis that the majority of the energy in so- 
lar wind turbulence has k±_ 3> feii (see IHowes" ]Eat|[2008a 
and references therein). In this case, there is a transi- 
tion to kinetic Alfven wave fluctuations at the scale of 
the ion Larmor radius. This Letter demonstrates that a 
dissipation range comprised of kinetic Alfven waves pro- 
duces a reduced fluctuating magnetic helicity signature 
consistent with observations. 
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using B 2 (fc) = B 2 {[(fc Pl )- 1/3 + (fcft) 4/3 ]/[l + (^IP- 
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